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ABSTRACT 

We present high sensitivity Very Long Baseline Interferometry (VLBI) obser- 
vations 806 days after the 7-ray burst of 2003 March 29 (GRB 030329). The 
angular diameter of the radio afterglow is measured to be 0.347 ± 0.09 mas, 
corresponding to 0.99 ± 0.26 pc at the redshift of GRB 030329 {z = 0.1685). 
The evolution of the image size favors a uniform external density over an R~ 2 
wind-like density profile (at distances of R > 10 18 cm from the source), although 
the latter cannot be ruled out yet. The current apparent expansion velocity of 
the image size is only mildly relativistic, suggesting a non-relativistic transition 
time of £nr ~ 1 yi. A rebrightening, or at least a significant flattening in the 
flux decay, is expected within the next several years as the counter-jet becomes 
visible (this has not yet been observed). An upper limit of < 1.9 c is set on the 
proper motion of the flux centroid. 

Subject headings: gamma rays: bursts 



1. Introduction 



Our understanding of the ori gin of gamma-ray bursts (G RBs) has advanced rapid ly in 
recent years since the fi rst X-ray (jCosta et a/.lll997l ). optical (Ivan Paradijs et a/.lll997l ) and 
radio (IFrail et al. 1119971 ) afterglows were discovered. We now know that a GRB represents 
the quick (< 1 2 s) release of a l arge amount of energy - about 10 51 erg, when corrected for 
jet collimation (IFrail et al.ll200ll ; lBloom. Frail fc Kulkarnill2003l ). The prompt 7-ray emission 
is usually attributed to internal shocks due to variability within the ultra-relativistic outflow 
(with a Lorentz factor > 10 2 ), while the afterglow is from a strong relativistic shock that is 
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driven into the external medium as it decelerates the ejecta (for a review see iPiranl 120051 ) . 
The afterg low emission is predominantly synchrotron radiatio n, and its spectral and temporal 
evolution ( Sari. Piran fc Narayan 1998 ; Granot fc Sari 20021 ) are governed by such factors a s 



the structure and dynamical evolution of the relativistic jet (for a review see Granotl 



2006 



and t he environment (e.g., constant density versus a wind-like density profile; I Chevalier fc Li 
20001 ) which are not well constrained. 



GRB 030329, discovered by HETE-2 OVanderspek et al.ll2003l ). and localized rapidly in 
the optical bands by Peterson & Price (20 03) represents a un ique opportunity for VLBI 
observations. At a redshift of z — 0.1685 (IGreiner et al.l 120031 ) this is the second nearest 
cosmo logical burst in t he northern sky detected to date (only GRB 060218 is nearer at 
a redshift z = 0.0331, iMirabal &; Halpernl 120061 ). Observation with the VLA show the 
GRB 030329 radio afterglow to be the brightest dete cted to date, with a maximum flux 
density of 55 mJy at 43 GHz one week after the burst (IBerger et al. 1120031 ) . 



Most VLBA studies so far have been restricted to rather distant GRBs (z ~ 1). While 
these observations have proved to be invaluable in studying the flux evolution of the afterglow 
and providing accurate p ositions (in one cas e localizing the GRB within 0.3 kpc of the center 
of its radio host galaxy; IBerger et all 120011 ). in all cases the afterglow faded before it could 
be resolved by the VLBA. In the case of GRB 030329 we have an exceptional instance of a 
radio afterglow near and bright enough that it can be resolved with VLBI. Its exceptional 
brightness enables the radio afterglow to be studied for much longer times than has been 
typical. 



The VLBA campaign on GRB 030329 (ITavlor et al\ 12004 120051 ) resolved the afterglow 
image and showed an increase in its diameter from 65 ± 22 /ias at t = 25 days after the burst 
to 172 ± 43 /ias at t = 83 days. Only slow growth was seen out to t = 217 days when the 
diameter was measured to be 176 ±80 /ias, possibly indicating a transition to non-relativistic 
expansion. Here we present very late time observations of GRB 030329 that were performed 
at t = 806 days after the burst with the goal of differentiating between afterglow models, and 
constraining the jet structure and dynamics as well as the environment. Such observations 
can also test the hypothesis of a double-sided jet, which would imply a significant contribution 
from the counter-jet as the flow becomes sub-relativistic at late times. 

Assuming a A cosmology with Hq = 71 km s -1 Mpc -1 , Qm = 0.27 and = 0.73, the 
angular-diameter distance of GRB 030329 at z = 0.1685 is d>A = 587 Mpc, and 1 milliarcsec 
corresponds to 2.85 pc. 
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2. Observations 

2.1. Late Time VLBI Observations 

Late time observations require high sensitivity to detect the fading afterglow. To achieve 
this high sensitivity we employed the largest aperture radio telescopes available, and used 
a bandwidth of 64 MHz in two circular polarizations. The VLBI observations were taken 
at 5 GHz on 2005 June 12, 806 days after the burst, with a global array including the 
100 m Green Bank Telescope (GBT) of the NRAO0, the Effelsberg 100 m telescopff, the 
305 m Arecibo telescopdfl, the Westerbork (WSRT) tied arrajfl, and the 25 m Markll tele- 
scope of Jodrell Banlo (replaced from the initially scheduled Lovell telescope due to the 
fast source-switching phase-reference schedule). In all, the combined effective collecting area 
is 0.097 km 2 . The on-source time on GRB 030329 was 1.6 hr for all telescopes except for 
Arecibo (0.9 hr). Data were recorded at 1 Gb s" 1 and 2-bit sampling at all stations except 
for Green Bank, where data were recorded at 512 Mb s _1 and 1-bit sampling. The corre- 
lator saw the Green Bank data as 2-bit data, with the magnitude bit always set to high. 
The post-correlation 2-bit van Vleck correction thereafter processed the data as usual. The 
observations were correlated at the Joint Institute for VLBI in Europe (JIVE). 

The nearby (0.8°) source J1048+2115 was used as the phase-reference calibrator, and 
was observed regularly in a cycle of 2 minutes on the target and 1 minute on the calibra- 
tor. Since J1048+2115 is not tied to the International Celestial Reference Frame (ICRF), 
the ICRF source J1051+2119 was observed every 30 minutes. This allowed an accurate 
determination of the position of the phase-reference source and subsequently the target. 
Self-calibration solutions for the amplitude were transferred from the stronger J1051+2119 
calibrator, which corrected for a large amplitude offset affecting all the left hand circular 
polarization IFs of the Markll telescope. In proposing these observations, we aimed for a 
theoretical noise of around 6 /iJy/beam. However, a number of factors resulted in a de- 
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creased sensitivity in our final map. These factors includes the replacement of Lovell with 
the smaller Markll telescope, and missing or bad data. As a result, a noise of 15 /zJy/beam 
is reached in the final image. 



2.2. VLA Observations 



Multi-frequency V LA observations of GRB 030329 have been presented ( IBerger et al. 
20031 ; iFrail et al\ 120051 ) between 1 day and 1 yr after the burst. Here we present VLA 



observations at 5 GHz taken from the NRAO archive between 380 and 740 days after the 
burst (see T able 1 for more det ails). The VLA data were reduced in AIPS following standard 
procedures (ITaylor et a/.lll9991 ). The nearby (1.6° separation) calibrator J1051+213 was used 
to determine the complex antenna gains. Absolute flux density calibration was tied to short 
observations of 3C286. 



Using this VLA data in combination with published WSRT data (Ivan der Horst et al. 



20051 ). we find that the fall in flux densities from day 59 to day 806 (Fig. 1) are well described 



by a power law, F v oc t~ a , with a slope a = 1.23 ± 0.03. In the end of 2009, the data 
recording rate of the VLBA is planned to increase to 4 Gb s _1 compared to its current 
limit of 0.5 Gb s™ 1 . The inclusion of the VLBA at this recording rate, as well as the large 
antennas used for this experiment, should provide an increase in sensitivity by a factor ~ 3, 
or a theoretical rms noise of ~ 2 /i Jy. Assuming that GRB 030329 continues to fade in the 
same way, then we expect a flux density in the end of 2009 (6.6 years after the burst) of 
~ 70 /iJy. Thus it should be quite practical to follow the growth of GRB 030329 for at least 
another 3-5 years. 



3. Results 

We fit a symmetric, two-dimensional circular Gaussian to the measured visibilities on 
G RB 030329 and find a full-width at half maximum (FWHM) size of 0.347 ± 0.09 mas. As 



in 



Taylor et al\ (120041 ) the error on the size is estimated from signal-to-noise ratios and from 
Monte-Carlo simulations of the data using identical (u,v) coverage, similar noise properties, 
and a Gaussian component of known size added. The standard deviation of the recovered 
sizes, model-fitted in the same way as we treat the observations, was found to be 0.07 mas. 
Given two similar estimates of the error, we adopt the more conservative limit of 0.09 mas 
based on signal-to-noise ratios and the resolution of the observations. 

We also obtain a position for GRB 030329 of R.A. 10 h 44 m 49.95957 s and Dec. 21°3L 17.4375" 
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with an uncertainty of 0.2 mas in each coordinate. 



4. Discussion 

4.1. Angular Size Measurements 

The entire history of expansion for GRB 030329 is shown in Fig. [2j The first measure- 
ment at 15 days com es from a model-dependent estimate of the quenching of the scintillation 



(IBerger et al. II2003I ). The uncertainties on this size estimate are large due to the dependence 
of the measurement on estimated properties of the interstellar medium along this particular 
line-of-sight. Furthermore, a limb-brightened image has a smaller effective area (similar to 
a thin ring) compared to a uniform image of the same diameter. This reduces the effective 
number of regions within the image with independent brightness variations for a given image 
diameter. Thus, in a limb-brightened image (as expected here), the image diameter might 
be under-estimated by a factor larger than that for the fit to a circular Gaussian that is 
used in the other epochs (see discussion below Eq. [I]), thus potentially causing a systematic 
under-estimate of the source size at 15 days, compared to the other epochs that are shown 
in Fig. H 

Our directly measured size of 0.347 ± 0.09 mas a t 806 days is a fac tor of ~ 2 larger 



than the size of 0.176 ±0.08 mas measured at 217 days (ITavlor et a/.ll2005f). w hich was quite 
close to the size of 0.172 ± 0.043 mas measured at 84 days (ITaylor et a/.ll2004l ). We consider 
the possibility of no expansion between days 217 and 806, which at a first look at Fig. [2] 
might be a plausibility. However, the points at 217 and 806 days are separated by 2a, and 
the edges of the la regions just touch each other. At the other extreme, the size could have 
changed with a factor 4.6. In this case, the slope will be inconsistent with the earlier point 
at day 84, unless the source increased its expansion speed or shrank in size. We find this 
implausible, and that an expansion equal to (or slower than) a factor of ~ 2 is more likely. 

Adopting measured size increase of a factor of ~ 2, this implies a continued growth of 
the image size, and confirms a gradual decrease in its apparent expansion speed, /3 app c. In 
Fig. 3 we show the decline of (/? app ) with time, calculated from the measured sizes using 

in V _ (j_+ Z )R± _ Qr<Im ,-,s 

UW- cf - ct > W 

where 6r = R±/d>A = (1 + z)R±/d.M and R± are the angular and physical radius of the 
image, respectively, and du is the proper distance to the source. We note that similar to 
arguments for no expansion in Fig. 2, it could be argued that the apparent expansion speed 
(Fig. 3) at later times is constant. This would mean that the source size grows linearly with 
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time (i.e. a factor of 3.7 growth between day 217 and 806). As mentioned above, this type 
of growth requires a rather strange behavior given the point at day 84. Thus, a decrease in 
the apparent expansion speed seems more reasonable. 

The values for the apparent size depend somewhat upon the intrinsic surface brightness 
profile, which for c ircular Gaussian, a uniform disk, and a rin g, vary roughly as 1:1.6:1.1 



( Taylor et q/.ll2004l). In practice, the image is limb-brightened (IGranot. Piran fc Saril 11999 



Granot fc Loebl l200ll ). and perhaps slightly closer to a thin ring than to a uniform disk, 
suggesting a correction factor of ~ 1.4. This correction factor might change with time, since 
the surface brightness profile within the afterglow image is expected to change across the non- 
relativistic transition, and possibly also during the post-jet break relativistic stage. However, 
the afterglow image during both of these stages has not been studied yet in sufficient detail 
to provide a reliable description of the exact temporal evolution of this correction factor. At 
any rate, this effect is expected to be fairly mild, and smaller than the factor of ~ 2 growth 
in the image size between 217 days and 806 days. 

Taking this correction factor into account, we find (/? apP ) ~ 1-2 after 806 days, i.e. 
that the current expansion velocity (which is expected to be somewhat smaller than the 
average value (/3 a pp)) is only mildly relativistic, while at earlier times it was considerably 
higher and more relativistic, as implied by the value of (/3 app ) ~ 8 after 25 days. The non- 
relativistic transition occurs when » pp pa 1, i.e. when (0 app ) ~ 2, the exa ct value depending 



on the detailed expansion history (IGranot. Ramirez- Ruiz fc Loebl 120051 ). This suggests a 
non-relativistic transit ion time of about £nr ~ 1 yr. Based on a joi nt temporal and spe ctral 
fit to the light curve, iFrail et all (120051 ) find £nr ~ 50 days, and iResmi et al\ (120051 ) use 
simultaneous radio, millimeter and optical observations to derive model-dependent non- 
relativistic transiti on times of £nr ~ 4 2 and ~ 63 days. Moreover, multi-frequency radio 
data presented by Ivan der Horst et al.l (120051 ) is fitted with a £nr ~ 80 days. All these 



transition times are condsiderably shorter than our estimate. In principle their models could 
be refined using the direct size estimates presented here to better constrain the energetics of 
the burst and other physical parameters of the models. 



Radio re-brightening of GRB 030329 has been predicted by IGranot &: Loebl (120031 ) as 
the count er jet becomes no n-relativistic and its emission is no longer strongly beamed away 
from us. iLi fc Song! (120041 ) estimated a level of 0.6 mJy at 15 GHz, 1.7 yr after the burst, 
due to this effect. The actual flux density around that time, however, was ~ 3.3 times lower 
than their prediction (0.36 mJy at 4.86 GHz after 621 days, with F u oc z/ -06 in the radio, 
between the two frequencies, implying pa 0.18 mJy at 15 GHz). More importantly, we find 
that the late time 4.9 GHz light curve, from 59 days and up to 806 days, is consistent with a 
single power law decay of F v oc t- L23±a03 ( se e Fig. [TJ, which implies that up to ~ 800 days 
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after GRB 030329 there is still no significant contribution to the observed flux from the 
counter jet. Light travel effects might cause the contribution from the counter jet to become 
significant at a somewhat later time, up to a factor of several after t N R when (3 apx , pa 1 , i.e. 



up to several years after the burst (this effect was taken into account by lLi & Songll2004 but 
they estimated a somewhat lower value for £nr of ~ 120 days). Alternatively, if the counter 
jet encounters a somewhat lower external density than the forward jet, it would become 
non-relativistic at a later lab frame time, and light travel effects would cause its contribution 
to the observed flux density to become significant at an even later time (as well as being 
intrinsically less luminous than the forward jet at the same radius). If such a re-brightening 
(or at least a significantly shallower flux decay) eventually occurs within the next several 
years, it may enable measurement of the image size for a longer time. 

Fig. H] demonstrates how the measured evolution of the image size for GRB 030329 
compares with the predictions of different theoretical models. The data constrain the external 
density profile at radii R > 10 18 cm, and favor a uniform density external medium (p ext = 
Ar~ k with k = 0) over a stellar wind environment (k = 2, as expected for a massive star 



progenitor; IChevalier fc Lill2000l ). although it is hard to rule out the latter due to the sizeable 
errors on the source size, and the theoretical uncertainty in the jet structure and dynamics. 
The degree of lateral spreading of the jet is harder to constrain, since its effect on the 
evolution of the observed image size around the times when it was measured is somewhat 
smaller than that of the external density profile (see Fig. H]). 



4.2. Constraints on Proper motions 

Solving for proper motion using all the VLBI observations to date, we derive fi T ^ = 
0.12 ± 0.15 mas yr _1 and /idee. = —0.11 ± 0.08 mas yr _1 , or an angular displacement over 
806 days of 0.36 ± 0.38 mas (Fig. 5). These observations are consistent with those reported 



by lTaylor et al\ (120041 120051 ) . and impose an even stronger limit on the proper motion. The 
implied limit on the proper motion in the plane of the sky is < 0.38 mas (corresponding to 
< 1.08 pc) in 806(1 + z)- 1 « 690 days, or < 1.87 c (2 a). 

Our new upper limits on the proper motion are about a factor of 2 better than pre- 
vious limits. The proper motion of the flux centroid is constrained to be less than the 
current diameter of the image, which is consistent with the expectations for a narrow dou- 
ble sided jetfl where the maximal angular displacement of the flux centroid is < o bs#NR 



6 The observed image angular diameter around £nr is expected to be > 2# -R(iNR) / 'd-A — 2#o$nr, which 
is a factor of > 29o/9 h s > 2 larger than the expected maximal angular displacement of the flux centroid. 
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and is obtained around £nr> where #nr = -R(^nr) / (IGranot fe Loebl 120031 ). since for 
GRB 030329 the viewing angle is within the jet, Q \^ < 9 . This consistency is indepen- 



(IGranot. Nakar & Piran 


2003; 


Bereer et al. 112003: 


Sheth et al. 


2003; 


Piran. Nakar & Granot 


2004; 


Granotl 


2005: 


Gorosabel et alhood). 



5. Conclusions 

We present measurements of the 5 GHz flux density and size of the radio afterglow 
from GRB 030329 over a period of 2.2 years. These observations clearly demonstrate that 
the expansion velocity has slowed down over time, with a transition to the non-relativistic 
regime at ~1 yr. The evolution of the afterglow image size favors a uniform external density 
over a wind-like stratified external medium (at R > 10 18 cm), although the data still do not 
enable to conclusively rule out a wind-like density profile (see Fig. @J. We also present a 
tight upper limit on the apparent motion of the flux centroid (< 1.87 c, a factor of 2 better 
than previous limits), which are consistent with the predictions of standard afterglow theory. 

Although GRB 030329 has faded considerably, it is likely that VLBI observations will 
continue to be useful for determining expansion size estimates in the future. With VLBI 
recording rates expected to increase to 4 Gb s _1 over the next few years, thermal noise levels 
of the VLB A at 5 GHz should decrease to ~ 8.5 /iJy, and arrays of larger antennas, such 
as the one used in this work, will reach theoretical thermal noise levels of ~ 2 /iJy. This 
anticipated increase in VLBI capability should keep pace with the flux density of GRB 030329 
for the next few years. A rebrightening, or at least a significantly shallower flux decay, is 
expected within the next several years, as the counter jet begins to significantly contribute 
to the observed flux. Further VLBI observations should improve constraints on the density 
profile of the external medium into which the afterglow shock is expanding, and should 
eventually reveal the shape of the afterglow. 

This research has made use of NASA's Astrophysics Data System. The research of JG 
is supported by the US Department of Energy under contract DEAC03-76SF00515. 
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Table 1. Observational Summary 



Date 


At 


Frequency 


Flux Density 


Instrument 




(days) 


(GHz) 


(A*Jy) 




2004 Apr. 13.11 


380.63 


4.86 


631 ± 50 


VLA-C 


2004 May. 29.12 


426.64 


4.86 


516 ± 39 


VLA-C 


2004 Jun. 01.15 


429.67 


4.86 


485 ± 44 


VLA-D 


2004 Jul. 19.94 


478.46 


4.86 


441 ± 43 


VLA-D 


2004 Aug. 15.79 


505.31 


4.86 


528 ± 53 


VLA-D 


2004 Sep. 03.79 


524.31 


4.86 


331 ± 28 


VLA-A 


2004 Dec. 09.66 


621.18 


4.86 


361 ± 36 


VLA-B 


2005 Jan. 29.38 


671.90 


4.86 


346 ± 51 


VLA-B 


2005 Apr. 07.21 


739.73 


4.86 


292 ± 37 


VLA-C 


2005 Jun. 12.78 


806.30 


4.99 


240 ± 17 


EB+GBT+WB+AR+JB 



*NOTE - EB = 100m Effelsberg telescope; GBT = 105m Green Bank Telescope; 
WB = phased Westerbork array; AR = 305m Arecibo telescope; JB = 25m Markll 
Telescope. 
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Fig. 1. — The 4.9 GHz light curve for GRB 030329 from publicly available data and these 
observations. The VLBI epochs are indicated by vertical lines between 3 days and 806 days 
after the burst. The data from day 59 to day 806 are well described by a power law, 
F v oc t-i-23±o.o3 ^ soM line y 
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Fig. 2. — The apparent expansion of GRB 030329 derived from measurements and limits 
on t he angular siz e as a function of time. The two upper limits at 7.7 days and 51 days are 
from lTaylor et al\ (120041 ) . as are the measurements o n days 25 and 83 (open squares). The 
measurement on day 217 is from iTaylor et al\ (120051 ) (open square). The measurement on 
da y 15 (filled cir c le) is a model dependent estimate based on the quenching of scintillation 
by Berger et al. (j2003 ). Finally, the measurement on day 806 (open square) comes from 
this work. The lower panel shows the evolution of the average apparent expansion velocity. 
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Fig. 3. — Evolution of the average apparent expansion velocity derived from direct size 
measurements, and assuming a Gaussian intrinsic surface brightness profile. 
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Fig. 4. — Tentative fits of theoreti cal models for the evolution of the source size (from 
Granot. Ramirez- Ruiz fc Loebll2005l ) to the observed image size (of diameter 2R±) of the 
radio afterglow of GRB 030329 up to 83 days, with the addition of the measured sizes at 
217 days (from Taylor et a/i I2OO5I ) . and at 806 days (from this work). In model 1 there is 
relativistic lateral spreading of the GRB jet in its local rest frame, while in model 2 there is 
no significant lateral expansion until the jet becomes non-relativistic. The external density 
is taken to be a power law with the distance r from the source, p cxt = Ar~ k , where k = 
for a uniform external density while k = 2 is expected for a stellar wind environment. 
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Fig. 5. — The positions derived from the observations in eight epochs relative to the first 
determination on April 1st at 8.4 GHz. Observations at multiple frequencies at a given 
epoch have been plotted separately since they are independent measurements. A circle with 
a radius of 0.26 mas (2 a) is shown to encompass all measurement s except those take n within 
the first eight days at 5 GHz, which suffer from systematic errors ( ITaylor et al\\2004 ). Taken 
together these observations provide a constraint on the proper motion of 0.35 ± 0.38 mas 
over 806 days. 



